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ABSTRACT 


A  parametric  study  of  a  body  of  revolution  is  conducted  utilizing  existing  semi- 
empirical  methods  for  the  calculation  of  hydrodynamic  coefficients.  The  geometry  of 
the  body  is  analyzed  in  non-dimensional  length  and  volume  parameters.  The  effects  of 
varying  the  nose,  mid-body,  and  base  fractions  of  the  body  on  the  hydrodynamic 
coefficients  are  generated  and  illustrated  graphically.  Equations  for  the  hydrodynamic 
coefficients  are  then  determined  from  the  non-dimensional  parameters.  The  results 
can  be  used  to  evaluate  fundamental  maneuvering  characteristics  early  in  the  design 
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I.  INTRODUCTION 


A.  GENERAL  CONCEPTS 

Ships  are  designed  to  operate  in  a  hostile  environment.  Even 
during  peacetime  the  ship  must  navigate  the  water  and  arrive  at  it’s 
destination  intact  and  ontime.  The  motion  of  the  ship  is  affected  by 
the  presence  of  other  ships,  hazards  to  navigation,  established  traffic 
patterns,  and  the  sea.  The  motion  of  the  ship  is  directed  and  the 
track  of  the  ship  is  plotted  on  the  chart.  The  response  of  the  ship  to 
changes  in  speed  and  rudder  commands  is  studied  under  standard 
assumptions  of  marine  vehicle  dynamics.  The  forces  and  moments 
subjected  to  the  ship  due  to  its  relative  motion  in  the  water  are  the 
basis  for  this  study.  Figure  1  illustrates  the  basic  concepts  and 
definitions. 

Y  Force  v 


Figure  1.  Overview  of  forces  and  velocities.  [1] 
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The  motion  of  the  ship  can  be  described  in  the  global  sense  and 
in  the  body  fixed  reference  frame.  The  body  fixed  reference  frame  is 
utilized  in  this  study.  The  ship's  velocity  along  its  x  axis  is  u.  The 
sideslip  velocity  is  v  and  the  angular  velocity  is  r. 

When  the  water  impacts  the  ship  at  an  angle  of  attack,  lift  and 
drag  forces  are  developed.  These  impart  a  lateral  Y  force  and  a 
pitching  N  moment.  We  can  determine  the  lateral  Y  force  resulting 
from  a  sideslip  velocity  v  and  an  angular  velocity  r.  We  can 
determine  the  pitching  N  moment  resulting  from  a  slideslip  velocity 
v  and  angular  velocity  r.  In  review  of  the  effects  of  acceleration, 

V  and  j-  ,  these  effects  can  be  computed. 

A  measure  of  the  ship's  stability  in  design  is  its  ability  to 
maintain  or  regain  straight  line  motion  when  forces  act  on  its  shape. 
The  linearized  models  for  equations  of  motion  utilize  hydrodynamic 
derivities  of  the  Y  force  and  the  N  moment  to  describe  their 
behavior.  There  are  a  number  of  mathematical  methods  for 
determining  the  hydrodynamic  derivatives  and  LT  Wolkerstofer  [2] 
summarized  them  in  his  thesis.  In  our  work,  the  semi-empirical 
methods  that  utilize  the  geometric  considerations  of  a  body  of 
revolution,  typical  of  a  modern  submarine  are  best  suited  for  our 
purpose. 

B.  PROGRAM  APPLICATION 

The  body  of  revolution  that  we  shall  consider  is  a  basic 
submarine  shape.  The  nose  is  elliptical,  the  mid-body  is  cylindrical, 
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and  the  base  is  conical.  This  is  an  approximate  SUBOFF  model  [2]  and 
it  is  also  the  shape  of  the  slice  pod  which  is  illustrated  in  Figure  2. 


Figure  2.  A  body  of  revolution. 

This  axisymmetric  shape  is  a  simple  model  but  the  concepts 
applied  to  this  model  can  be  broadened  to  fit  any  shape.  Our  goal  is 
to  take  a  shape  that  has  been  tested  and  to  modify  its  parameters  to 
see  how  geometric  changes  effect  the  hydrodynamic  derivatives 
which  affect  the  stability  and  turning  characteristics  of  the  platform. 
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H.  DETERMINATION  OF  HYDRODYNAMIC  COEFFICIENTS 


A.  THE  GENERAL  PHYSICAL  DESCRIPTION 

The  body  of  revolution  consists  of  three  main  sections.  The 
forward  section  is  an  ellipsoid,  the  mid-body  section  is  a  cylinder  and 
the  base  section  is  a  cone.  This  is  the  shape  utilized  for  numerous 
studies  and  is  the  shape  of  the  slice  pod.  The  generic  figure  is 
illustrated  in  Figure  3. 

Nose:  Ellipise  Parallel  Mid-body:  Cylinder  Base:  cone  shape 


#:  Individual  section  geometric  center 
Figure  3.  Geometric  parameter  description  for  body  of  revolution. 

The  various  methods  utilized  for  the  estimation  of 
hydodynamic  coefficients  were  compared  by  LT  Wolkerstofer  [2]. 
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The  body  of  revolution  hydrodynamic  coefficients  depend  on 
semi-empirical  relations  to  account  for  the  viscous  or  vortex  effects 
on  the  body.  The  methods  based  primarily  on  geometric 
considerations  are  applied  in  this  parametric  study.  The 
hydrodynamic  coefficients  will  be  determined  from  the  United  States 
Air  Force  Data  Compendium  (USAF  DATCOM)  method  which  is 
summarized  by  Peterson  [3].  The  acceleration  hydrodynamic 
coefficients  will  be  determined  by  the  methods  of  Humpreys  and 
Watkinson  [4]. 

To  initiate  the  parametric  study  we  needed  to 
non-dimensionalize  a  number  of  characteristics  of  the  body  of 
revolution.  The  lengths  of  each  section  of  the  body  are  described  by 
the  fractional  amount  of  the  total  length  as  in  the  nose  fraction. 


The  slenderness  ratio  is  defined  using  the  total  length. 

„  h 

S  =  -j  (2) 

a 

The  volume  the  body  of  revolution  would  be  determined  from 
the  sum  of  the  three  sections  using  standard  volume  formulas  but 
with  all  the  lengths  referenced  to  the  total  length  of  the  body  via 
fractional  amounts. 

v={  ^](2F»+3f»+f*)  <3> 

V  1Z 


6 


From  this  relation  it  is  clear  that  the  mid-body  fraction  is  the 
dominant  term  in  determination  of  the  volume  amount.  This  is 
consistent  with  the  realization  that  the  mid-body  diameter  is  the 
maximum  diameter  for  the  body  and  constant  mid-body  length.  The 
determination  of  the  volume  is  important  because  it  is  the  volume 
which  determines  the  buoyant  forces  resulting  from  the  shape.  A 
larger  volume  will  have  a  larger  buoyant  force. 

When  the  volume  is  non-dimensionalized  the  diameter  of  the 
body  of  revolution  will  be  determined  in  the  following  relation. 

d  =  V(12V)/K)(2f-.+3f»+fs)  <4> 

The  geometric  center  of  the  body  of  revolution  will  depend  on 
the  lengths  of  the  individual  sections.  We  utilized  standard  formulas 
for  each  section  to  find  the  geometric  center  for  each  section  and 
then  combined  then  to  determine  the  offset  from  the  geometric 
midpoint  of  the  mid-body  section. 


This  offset  from  the  center  of  the  mid-body  section  is  then  applied  to 
find  the  geometric  center  for  the  body  which  is  the  same  as  the 
center  of  bouyancy  for  the  body  in  horizontal  motion. 

L  =  L+If~xb  <«) 
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The  axial  position  where  the  flow  becomes  predominantly 
viscous  is  a  function  of  the  overall  length  determined  from  the  point 
on  the  body  of  maximum  slope  [3].  The  axial  position,  Ly  is  utilized 
for  the  determination  of  a  few  parameters  affecting  the 
hydrodynamic  coefficients. 

lv  =  0.905  xZB  (7) 

The  geometry  of  the  body  at  Lv  is  illustrated  in  Figure  4. 


Sv:  Cross-section 
Area  of  Base 


Stb:  Profile  Area  of  Base 


LB-Lv 


Figure  4.  Sy  and  Stb  illustration. 


The  maximum  profile  area  is  determined  at  the  mid-body 
section  where  the  diameter  is  maximum. 


(8) 


The  radius  of  the  body  at  the  axial  position  Ly  is  determined 
by  realizing  the  base  profile  is  a  triangle. 


The  cross  sectional  area  at  the  axial  position  is  therefore 


Sv  =  nr2  (10) 

and  the  profile  area  of  the  base  at  the  axial  position  is  as  follows: 

Slb=r(lB-lc)  (11) 

The  drag  coefficient  for  the  body  is  initially  established  by 
Fidler  and  Smith  [5]  and  utilized  by  Wolkerstofer  [2]-  The  value  of 
the  drag  coefficient  will  be  modified  within  the  parametric  study  in 
Chapter  III. 

Cd  =  0.29  (12) 

B.  HYDRODYNAMIC  COEFFICIENT  CALCULATION 

From  the  basic  physical  description  of  the  system  the 
hydrodynamic  parameters  can  now  be  determined. 

The  Lamb's  coefficients  (ki  and  k2)  of  inertia  are  for  a  prolate 
ellipsoid  in  axial  and  cross  flow.  These  parameters  affect  the 
lift/angle  of  attack  curve  slope  [3]. 


C„  =  2  s. 


k2  -  k1 


(13) 


The  hydrodynamic  coefficient  Yv  is  the  normal  force/angle  of 
attack  curve  slope  [3]. 


vprtme  1 2 

V 


(C„+Cd) 


(14) 


The  pitching  moment/angle  of  attack  curve  slope  was 
simplified  using  MATHCAD  and  its  symbolic  manipulator  and  is 
based  on  the  cross  sectional  area  of  the  body  of  revolution  [3]. 


_  2(A:2  kx ) 


-ra*2(3 -/„)  n rf2 
12  ml 

r-3xJ2v+6lvxJb  +  2l3v  -  31  fl2v-3lbl2v 
^+6xmlvlf  ~6xmljlb  +3 ljlb  -3 xmlf  +1^ 


(15) 


The  hydrodynamic  coefficient  Ny  is  the  pitching  moment 
coefficient  [3]. 


N  -C  x-^ 

vprime  ma  j2 

lB 


(16) 


to 


The  lift/pitch  rate  curve  slope  is  a  function  of  the  lift/angle  of 
attack  curve  slope  [3]. 


r  —  r  1  _  Zul 

lq  -  '-la  1  * 

V  lB 


(17) 


The  rotary  derivative  hydrodynamic  coefficient  Yr  is  the 
normal  force/pitch  rate  coefficient  [3], 


c„sb 


(18) 


The  pitching  moment/pitch  rate  curve  slope  [3], 


X 

2  —  m 


V(lcb-xm) 


1_^L 

L 


S  l2 

°tblB 


Sfi bh 


(19) 


The  rotary  derivative  hydrodynamic  coefficient  Nr  is  pitching 
moment/pitch  rate  coefficient  [3]. 


CmqSb 


(20) 


The  acceleration  hydrodynamic  coefficients  Yy'  and  Ny'  are 
based  on  the  work  of  Humphrys  and  Watkinson  [4]. 
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We  computed  the  values  of  the  hydrodynamic  coefficients 
using  a  MATLAB  program  and  utilized  the  generalizations  that  will 
be  required  for  the  parametric  study.  The  results  were  compared  to 
the  DATCOM  SUBOFF  data  collected  by  LT  Wolkerstofer  [2]  and  the 
results  are  listed  in  Table  1. 
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Hydrodynamic 

Coefficient 

DATCOM 

Program 

Function 

Percent 

Error 

Yv' 

-0.0058 

-0.0058 

0 

Nv' 

-0.0136 

-0.0136 

0 

Yr' 

-0.0014 

-0.0014 

0 

Nr’ 

-0.0012 

-0.0011 

8.3 

YV’ 

-0.0153 

-0.0152 

0.6 

Nv’ 

0.0153 

0.0152 

0.6 

Nr' 

-0.0007 

-0.0007 

0 

Table  1.  DATCOM  and  program  comparison. 

The  results  of  the  program  agree  with  the  data  compiled  from 
the  actual  model  tested  using  the  semi-empirical  methods.  The 
modifications  to  simplify  the  analysis  based  on  the  given  shapes 
provided  accurate  results.  The  magnitude  of  the  Nr'  error  is  due  to 
the  small  value  of  the  Nr'  coefficient. 
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III.  PARAMETRIC  STUDY 


A.  INTRODUCTION 

In  conducting  a  parametric  study  of  the  body  of  revolution,  we 
needed  to  decide  which  parameters  to  vary,  which  to  hold  constant 
and  what  effects  these  variations  would  have  on  the  model.  The 
primary  goal  was  to  find  non-dimensional  parameters  that  would 
define  the  shape  sufficiently  to  accurately  predict  the  hydrodynamic 
coefficients.  The  non-dimensional  parameters  would  therefore  not 
be  constrained  to  a  given  particular  size  and  shape. 

The  data  which  completely  describes  the  body  of  revolution 
include  the  sectional  lengths  and  nominal  (mid-body)  diameter. 

From  this  the  volume,  cross  sectional  areas,  and  geometric  centers 
could  be  calculated  and  the  resulting  hydrodynamic  coefficients. 

This  initial  program  was  limited  to  specific  non-variable  inputs  only 
which  produced  single  case  results. 

The  initial  steps  in  broadening  the  program  were  due  to 
non-dimensionalizing  the  sectional  lengths  while  still  specifying  the 
overall  length  of  the  body.  The  diameter  was  determined  from  the 
slenderness  ratio.  In  describing  the  body  in  this  manner,  we 
determined  that  the  hydrodynamic  coefficients  were  independent  of 
the  specific  overall  length  of  the  body  as  the  bodies  parameters  were 
all  proportional  to  the  length.  We  could  therefore,  alter  the  sectional 
fractional  values  to  determine  the  effects  on  the  hydrodynamic 
coefficients. 
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In  order  to  study  these  variations  in  a  controlled  manner  we 
decided  to  hold  the  volume  of  the  body  as  a  constant  initially.  With 
the  volume  as  a  constant,  describing  the  sectional  fractional  values 
would  determine  the  remaining  parameters.  In  addition,  only  the 
nose  and  mid-body  fractions  would  be  altered  since  the  base  fraction 
would  also  be  defined  by  default. 

From  this  basis  we  pursued  two  approaches,  one  involving  a 
constant  diameter  and  varying  lengths  and  the  second  involving  a 
constant  length  and  varying  diameter.  This  vectoral  approach  was 
later  modified  to  allow  varying  lengths  utilizing  the  non-dimensional 
volume/length^  parameter. 

B.  PHYSICAL  REVIEW  OF  PARAMETER  MODIFICATIONS 

In  order  to  appreciate  the  multitude  of  possible  parameter 
modifications  Figures  5  through  8  are  presented  to  solidify  the 
physical  changes  incorporated  in  the  parametric  study.  Figure  5  is 
the  basic  shape  of  the  body  of  revolution.  We  have  repeated  this 
shape  behind  the  modified  shapes  for  Figures  6  through  8. 


Figure  5.  Basic  shape. 
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Figure  6.  Lower  Fn  and  slenderness  ratio  (shorter  and  fatter). 


Figure  7.  Higher  Fm  and  lower  Fb,  higher  slenderness  ratio(slimmer). 


Figure  8.  Lower  Fn,  higher  Fm,  lower  Fb,  constant  slenderness  ratio. 

C.  COMPUTER  MODEL  DEVELOPMENT 

In  the  process  of  the  parametric  study,  the  geometry  of  the 
body  would  be  altered  which  would  affect  the  assumptions  and 
variables  of  the  semi-empirical  method  utilized.  The  primary 
concern  was  with  the  drag  coefficient.  The  changed  shape  will 
induce  a  different  drag.  We  researched  the  significance  and 


magnitude  of  these  modifications  on  the  hydrodynamic  coefficients. 
The  value  of  drag  coefficient  was  a  constant  0.29  for  various  nose 
and  base  configurations  in  Fidler  and  Smith  [5]  but  at  the  extreme 
sectional  fractions  it  must  be  modified.  We  utilized  the  values  from 
Herner  [6]  in  the  extreme  cases.  When  the  nose  length  to  diameter 
ratio  was  less  than  1  the  coefficient  of  drag  was  smoothly  increased 
to  the  maximum  value  of  0.82.  When  the  base  length  to  diameter 
ratio  was  less  than  1  the  coefficient  of  drag  was  smoothly  increased 
to  the  maximum  value  of  0.64.  When  incorporated  into  the  program 
this  parameter  variation  had  minimal  effect  on  the  final  result  since 
the  other  components  prescribed  by  the  parameter  values  dominated 
the  solution  trend  and  our  specific  area  of  interest  is  not  at  the 
sectional  fractional  extreme. 

The  data  format  illustrated  in  Figure  9  demonstrates  the  limits 
of  the  model  and  is  useful  in  interpreting  the  mesh  graphs. 

Since  the  nose,  mid-body  and  base  fractions  together  comprise 
100%  of  the  body  the  values  in  the  lower  right  hand  side  diagonal 
are  not  possible  combinations  and  an  arbitrary  constant  value  is 
utilized  to  complete  the  matrix  and  provide  a  visible  floor  to  the 
mesh  graphs.  Appendix  A  is  the  MATLAB  program  that  generated 
the  matrix  data  with  the  full  range  of  nose  and  mid-body  fractons. 
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Mid-Body  Fraction 
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Matrix  Body:  Base  Fraction 
Figure  9.  Matrix  interpretation  for  data  analysis. 


D.  HYDRODYNAMIC  COEFFICIENT  MESH  GRAPHS 


The  Figures  10  through  17  were  generated  with  an  automatic 
scaling  feature.  This  allows  for  a  rough  comparison  of  the  magnitude 
of  each  parameter  to  be  compared  to  another  parameter.  The  Yv' 
variations  are  significant  while  the  acceleration  hydrodynamic 
coefficients  variations  in  magnitude  are  relatively  small.  The  graphs 
of  Ny'  and  Nr'  included  a  clipping  feature  to  allow  for  the  three 
dimensional  viewing.  The  clipping  only  involves  truncating  the 
values  after  the  solution  begins  to  reach  an  extreme  value. 

The  Figures  19  through  32  are  scaled  to  determine  the  shape  of 
the  hydrodynamic  coefficient  and  its  variations  over  the  entire  range 
of  sectional  fractions.  Each  hydrodynamic  coefficient  mesh  graph  is 
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coupled  with  a  two  dimensional  graph  which  is  not  clipped  and 
illustrates  the  varying  parameter  in  a  simpler  and  therefore  more 
clean  environment. 
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Figure  10.  Slenderness  ratio  mesh  graph. 
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Figure  11.  Yvprime  standard  scale  mesh  graph. 
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Figure  12.  Nvprime  standard  scale  mesh  graph. 
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Figure  13.  Yrprime  standard  scale  mesh  graph. 
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Figure  14.  Nrprime  standard  scale  mesh  graph. 
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Figure  15.  Yvdprime  standard  scale  mesh  graph. 
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Figure  17.  Nrdprime  standard  scale  mesh  graph. 


28 


Yvprime 


Nose  Fraction 


Mid-body  Fraction 


Figure  18.  Yvprime  close-up  scale  mesh  graph. 
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Figure  20.  Nvprime  close-up  scale  mesh  graph. 


Figure  21.  Nvprime  variations  with  mid-body  fraction. 


Figure  24.  Nrprime  close-up  scale  mesh  graph. 
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Figure  25.  Nrprime  variations  with  mid-body  fraction. 
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Figure  28.  Nvdprime  close-up  scale  mesh  graph. 
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Figure  29.  Nvdprime  variations  with  mid-body  fraction. 


E.  VOLUME  AND  LENGTH  VARIATIONS 


The  mesh  graphs  of  the  hydrodynamic  coefficients  were 
generated  holding  the  volume  of  the  body  of  revolution  constant.  We 
now  explored  variations  in  the  volume  and  length  parameters  of  the 
body  of  revolution  to  see  the  effects  on  the  hydrodynamic 
coefficients  as  illustrated  in  Figure  32. 


Figure  32.  Variations  of  V/L3  ratio  drawing. 

We  took  the  basic  body  and  coupled  the  volume  and  length  to 
determine  if  a  non-dimensional  volume/length3  ratio  was  valid.  The 
non-dimensional  results  matched  the  original  results.  The  value  of 
volume/length3  for  the  body  was  8.023xl0-3.  With  this  value  as  an 
average  value  we  varied  the  volume/length3  ratio  from 
6  to  lOxlO-3.  The  range  is  noted  on  all  figures  as  6  to  10  for 
simplicity.  The  smaller  value  is  due  to  a  smaller  volume  or  a  larger 
overall  length.  Figures  33  through  46  investigate  these  variations  at 
nose  fractions  of  10%  and  20%.  For  the  vast  majority  of  the  entire 
range  variations  of  volume/length3  result  in  linear  trends  in  all 
cases. 
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Figure  35.  Nvprime  V/L^  variations  at  10%  nose  fraction. 


Figure  36.  Nvprime  V/L^  variations  at  20%  nose  fraction. 
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V/l3  variations  at  20%  nose  fraction. 
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Figure  39.  Nrprime  V/L^  variations  at  10%  nose  fraction. 
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Figure  40.  Nrprime  V/L^  variations  at  20%  nose  fraction. 
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3.  Nvdprime  V/L^  variations  at  10%  nose  fraction. 


Figure  44.  Nvdprime  V/L^  variations  at  20%  nose  fraction. 
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Figure  45.  Nrdprime  V/L^  variations  at  10%  nose  fraction. 
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Figure  46.  Nrdprime  V/L^  variations  at  20%  nose  fraction. 
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IV.  DETERMINATION  OF  FUNCTIONAL  COEFFICIENTS 


A.  INTRODUCTION 


The  hydrodynamic  coefficients  for  various  combinations  of 
sectional  fractions  and  volume/length  ratios  have  now  been 
computed  and  graphed.  The  surface  functions  were  generated  using 
the  semi-empirical  methods  under  the  following  general  formula. 


HC  =  F 


V\ 

F  F  — 

1  n'1  m'  r 3 

L 


(25) 


The  goal  in  predicting  the  hydrodynamic  coefficients  would  be 
achieved  if  we  could  determine  the  relationship  along  the  surface 
function  based  on  the  variable  parameters. 

In  order  to  minimize  error  and  simplify  the  functional 
relationships  we  limited  the  range  of  the  variable  parameters  and  in 
each  case  the  range  of  the  parameters  covers  the  majority  of  the 
cases  of  interest.  The  nose  fraction  range  is  5%  to  25%.  The 
mid-body  fraction  is  40%  to  60%.  The  volume/length^  range  is 
6  to  10xl0‘3. 

In  reviewing  Figures  18-31,  in  our  range  of  interest,  the  curves 
appear  to  correspond  to  a  2n^  degree  polynomial.  In  reviewing 
Figures  33-46,  in  our  range  of  interest,  the  curves  have  a  definite 
linear  relationship.  The  difficulty  would  be  in  determining  the 
surface  functional  relationship  based  on  two  parameters  since 
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currently  there  is  no  built-in  subroutine  in  the  MATLAB  program 
toolbox  in  use  to  solve  this  problem. 

B.  PREDICTION  EQUATION 


Hard  copy  updates  and  distribution  of  improvements  in 
software  always  involves  time  delays.  In  corresponding  with  the 
Mathworks  company  I  learned  that  a  collection  of  m-files  existed  on 
the  internet.  I  utilized  some  of  these  m-files  in  the  development  of 
my  graphs  (the  legend  command  is  not  included  on  the  Macintosh 
MATLAB  professional  version).  I  downloaded  a  surface  equation 
least  squares  curve  fitting  program  that  I  could  modify  to  determine 
the  coefficients  of  the  surface  function.  I  tested  the  program  out 
using  a  number  of  different  surface  functions  and  concluded  the 
program  accurately  predicted  the  correct  polynomial  function. 
Appendix  B  includes  a  sample  program. 

With  the  functional  coefficients  determined  the  hydrodynamic 
coefficient  prediction  equation  would  be  as  follows: 


HC  =  AjF2n  +  A2FnFm  +  A3F2m  +  A4Fn  +  AsFm  +  A6 

+A7(l!~Cl) 


(26) 


The  coefficients  for  the  equations  are  listed  in  Table  2  for  the 
nose,  mid-body  fraction  and  volume/length  3  ratio.  The  constant  Cl 
is  8. 023x1  O' 3  and  is  the  nominal  value  for  volume/length^  ratio. 
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Figures  47  through  60  compare  the  theoretical  surface 
functions  to  the  predicted  surface  functions  using  the  same  axial 
scaling.  Figures  60  through  67  are  included  to  demonstrate  the 
percentage  error  of  the  prediction  equation  from  the  theoretical 
surface  function.  In  general,  the  percentage  error  is  small  with  the 
exception  of  the  Yr'  percentage  error  but  that  is  due  to  the  relatively 
smaller  magnitude  of  the  Yr  hydrodynamic  coefficient.  The 
acceleration  hydrodynamic  prediction  equations  are  very  accurate. 


H2  A 1 

A2 

A3 

A4 

A5 

A6 

A7 

Yy’  -0.0641 

-0.0641 

-0.0632 

0.0670 

0.0732 

-0.0263 

-0.5769 

Ny'  0.0277 

0.0499 

0.0266 

-0.0283 

-0.0301 

-0.0056 

-1.6357 

ram 

-0.0559 

-0.0292 

0.0310 

0.3160 

-0.0091 

-0.0880 

Nr  -0.0003 

0.0040 

0.0027 

-0.0012 

-0.0045 

0.0006 

-0.1590 

Yv'  0.0002 

0.0007 

0.0007 

-0.0008 

-0.0016 

-0.0144 

-1.8067 

Nv'  -0.0002 

-0.0007 

-0.0007 

0.0008 

0.0016 

0.0144 

1.8067 

Nr'  -0.0031 

-0.0046 

-0.0021 

0.0031 

0.0024 

-0.0013 

-0.0808 

Table  2.  Functional  Coefficients  for  prediction  equation. 

The  coefficients  of  Table  2  are  calculated  for  the  following 
nominal  values:  nose  fraction  15%,  mid-body  fraction  50%,  and  the 
volume/length  3  ratio  8.023x10*3.  These  parameters  represent  the 
mid-range  values. 
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Figure  47.  Yvprime  theoretical  mesh  graph. 
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Figure  48.  Yvprime  predicted  mesh  graph. 
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Figure  49.  Nvprime  theoretical  mesh  graph. 
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Figure  50.  Nvprime  predicted  mesh  graph. 
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Nrprime  Calculation 


Figure  53.  Nrprime  theoretical  mesh  graph. 


Figure  54.  Nrprime  predicted  mesh  graph 
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Figure  59.  Nrdprime  theoretical  mesh  graph. 
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Figure  60.  Nrdprime  predicted  mesh  graph. 
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Figure  61.  Yvprime  equation  percentage  error. 
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Figure  62.  Nvprime  equation  percentage  error 


Nrprime  Percent  Error 


Nose  Fraction  Mid-body  Fraction 

Figure  64.  Nrprime  equation  percentage  error. 
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V.  CONCLUSIONS  AND  RECOMMENDATIONS 


A.  CONCLUSIONS 

1.  Summary  Review  of  Parameter  Variations 

The  matrix  variations  of  the  nose  and  mid-body  fractions 
resulted  in  a  mixture  of  surface  functions  for  the  hydrodynamic 
coefficients.  There  was  not  an  overriding  trend  that  was  consistent 
with  all  the  hydrodynamic  coefficients.  The  trends  were  all  2nd 
order  polynomials  for  sectional  fractions  and  linear  for 
volume/length  ratio  variations.  The  larger  prediction  error  for  the 
Yr'  is  due  to  the  relatively  smaller  magnitude  variation  in  the 
hydrodynamic  coefficient.  The  acceleration  hydrodynamic 
coefficients  were  well  behaved  functions  and  the  prediction 
equations  are  highly  accurate.  The  acceleration  hydrodynamic 
coefficients  also  did  not  vary  greatly  in  magnitude.  Therefore,  the 
designer  can  alter  the  sectional  fractions  without  too  great  a  concern 
for  the  effects  due  to  acceleration. 

A  lower  volume/length  value  reduces  the  magnitude  of  the 
hydrodynamic  coefficient  and  this  is  consistent  for  all  the  cases 
studied.  For  a  constant  volume,  a  lower  volume/length  implies  a 
longer  slimmer  body,  and  a  higher  slenderness  ratio.  This  would 
imply  the  effect  of  the  diameter  on  the  body  has  a  greater  impact 
than  the  length. 
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2. 


Final  Comparison 


We  compared  the  program  results  to  the  expected  parameters 
of  the  SUBOFF  Body  and  the  Slice  Pod  The  results  are  listed  in 
Table  3. 


SUBOFF  Body _ Slice  Pod 


HC 

Predicted 

%  Error 

Predicted 

%  Error 

Yv' 

-0.0058 

-0.0054 

6.8 

-0.0141 

-0.0153 

7.8 

■Ml 

-0.0136 

-0.0140 

2.8 

-0.0395 

-0.0410 

3.6 

Yr' 

-0.0014 

-0.0008 

42.8 

-0.0014 

-0.0021 

30.0 

Nr’ 

-0.0015 

-0.0018 

18.2 

-0.0035 

-0.0034 

2.9 

Yv’ 

-0.0152 

-0.0151 

0.7 

-0.0442 

-0.0458 

3.5 

■SHM 

0.0152 

0.0151 

0.7 

0.0442 

0.0458 

3.5 

-0.0007 

0.0006 

14.2 

-0.0018 

-0.0019 

5.2 

Table  3.  Program  Comparison. 

Overall,  the  comparison  is  acceptable  with  the  understanding 
that  the  percentage  error  with  Nr'  is  due  to  its  relatively  small  value. 
The  Yr  percentage  error  is  within  the  tolerance  of  the  program. 


B.  RECOMMENDATIONS 


The  parametric  study  of  the  body  of  revolution  while  singular 
in  scope  presents  the  ability  to  pursue  variations  in  the  design  of 
components  and  their  effect  on  the  hydrodynamic  coefficients. 
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Recommendations  for  further  research  in  this  area  are  as  follows: 


•  Modify  the  program  to  evaluate  different  shapes. 

•  Evaluate  motion  in  the  vertical  plane. 

•  Calculate  stability  criteria  and  general  maneuvering 
performance  for  the  surface  mesh  functions. 
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APPENDIX  A.  MATLAB  ROUTINE  FOR  DETERMINING 
HYDRODYNAMIC  COEFFICINETS  WITH  VARYING  NOSE 
AND  MID-BODY  FRACTIONS 


A.  VDCURVES.M  PROGRAM 

%  LCDR  Eric  Holmes 
%  Constant:  Volume/LengthA3  ratio 
%  Thesis 

clear 

clear  global 

%  Basic  Parameters 

%  IB  =  Total  body  length 
%  In  =  Nose  length 
%  lm  =  Middle  body  length 
%  lb  =  Base  body  length 
%  d  =  diameter 
%  If  =  forebody  length  (lb  +  lm) 

%  Fn  =  Nose  fractional  length  of  total  body  length 
%  Fm  =  Mid-body  fractional  length  of  total  body  length 
%  r_Fm  =  maximum  fraction  of  mid-body 
%  Fb  =  Base  fractional  length  of  total  body  length 
%  N,I,k,t  =  index  counters 

%  Sb  =  maximum  cross  sectional  area  of  body  (assumimg  diameter 
%  d  and  length  IB) 

%  V  =  Volume  of  the  body 
%  xm  =  Geometric  middle  of  body 

%  xcb  =  Geometric  offset  from  the  center  of  the  middle  body 

%  for  gravity 

%  lcb  =  Geometric  center  of  gravity 

%  rho  =  density  of  water 

%  e  =  Munk  coefficient 

%  Bo  =  Munk  coefficient 

%  Ao  =  Munk  coefficient 

%  kl  =  Lamb's  inertial  coefficient 

%  k2  =  Lamb's  inertial  coefficient 

%  kb  =  Lamb's  inertial  coefficient 

%  Cdo  =  Drag  coefficient  at  zero  angle  of  attack 
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%  lv  =  length  where  viscous  flow  dominates 

%  r  =  radious  at  length  lv 

%  R  =  range  for  graphs 

%  S  =  Slenderness  ratio  matrix 

%  Sratio  =  slenderness  ratio 

%  Sv  =  cross  sectional  area  at  length  lv 

%  Stb  =  platform  sectional  area  at  lenght  lv 

%  Cla  =  Lift/Angle  of  attack  curve  slope 

%  Cma  =  Pitching  moment/ Angle  of  attack  curve  slope 

%  Cmq  =  Pitching  moment/pitch  rate  curve  slope 

%  Clq  =  Lift/Pitch  rate  curve  slope 

%  Iydf  =  mass  moment  of  inertia  of  displaced  fluid  mass 
%  Iydfjn  =  nose  body  component  of  Iydf 
%  Iydf_lm  =  middle  body  component  of  Iydf 
%  Iydf_lb  =  base  body  component  of  Iydf 

%  Hydrodynamic  coefficients 
%  Yvprime  =  normal  force  coefficient 
%  Nvprime  =  pitching  moment  coefficient 
%  Yrprime  =  normal  force/pitch  rate  coefficient 
%  Nrprime  =  pitching  moment/pitch  rate  coefficient 
%  Zwdotprime  =  acceleration  coefficient  (axisymetric  to  Yvdotprime) 
%  Yvdotprime  =  acceleration  coefficient  along  y  axis 
%  Nvdotprime  =  acceleration  coefficient  causing  yawing  moments 
%  Nrdotprime  =  acceleration  coefficient  in  addition  to  Iz 

%  Yvp,Nvp,Yrp,Nrp,Yvdp,Nvdp,Nrdp  are  hydrodynamic 
%  coefficient  matrices  each  combining  the  values  of  their  respective 
%  coefficient:  Yvpprime  =>  Yvp 

%  Yvp_d,Nvp_d,Yrp_d,Nrp_d,Yvdp_dNvp_d,Nrp_d  are  hydrodynamic 
%  coefficient  matrices  maintained  as  raw  data. 

%  (The  matrix  filler  is  -1.0) 

%  Yvp_s,Nvp_s,Yrp_s,Nrp_s,Yvdp_s,Nvp_s,Nrp_s  are  hydrodynamic 
%  coefficient  matrices  that  have  been  clipped  for  graphical 
%  presentation 

%  Initializing  data  and  empty  matrices 


S  =  [];  Yvp  =  [];  Nvp  =  [];  Yrp  =  [];  Nrp  =  [];  Yvdp  =[];  Nvdp  =  [];  Nrdp  =  []; 


R  =  (0:0.01:1); 
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v  =  1; 

IB  =  4.9941; 
rho  =  62.4/32.174; 

%  Drag  Function  Development 

ln_d  =  [0  .5  1  1.5  2  2.5  3]; 

cdjnd  =  [0.82  .45  .29  .29  .29  .29  .29]; 

cln  =  polyfit(ln_d,cd_lnd,4); 

lb_d  =  [0  .5  1  1.5  2  2.5  3]; 

cdjbd  =  [0.64  .38  .29  .29  .29  .29  .29]; 

clb  =  polyfit(lb_d,cd_lbd,4); 


%  Establishing  fractional  values  for  the  body 

Fn  =  (0:0.01:1); 
for  k  =  1:101; 
r_Fm  =  1  -  Fn(k); 

Fm  =  (0:0.01  :r_Fm); 

Fb  =  1  -  Fn(k)  -  Fm; 

%  Diameter  calculation 

d  =  ((12*V)./((pi*lB).*(2*Fn(k)  +  3*Fm  +  Fb))).A0.5; 

%  Length  calculation 

In  =  Fn(k).*lB; 
lm  =  Fm.*lB; 
lb  =  Fb.*lB; 

If  =  In  +  lm; 

%  Slenderness  ratio  calculation 

Sratio  =  lB./d; 

S  =  [S  Sratio  zeros(l,k-l)]; 

%  Body  calculation 
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Sb  =  (pi.*d.A2)/4; 
xm  =  IB/2; 


global  d  IB  If  In  lm  lb  xm  rho  I 
%  Center  of  gravity  and  bouyancy 

xcb  =  (pi*d.A2./12).*(2.*ln.*(lm/2  +  3*ln/8)  -  lb.*(lm/2  +  3*lb/4))/V; 

lcb  =  In  +  lm/2  -  xcb; 

%  Munk  coefficients 

e  =  (2./lB).*(((lB.A2)/4  -  Sb/pi))  A0.5; 

Bo  =  (l./e  A2)  -  ((1  -  e.A2)./(2*e.A3)).*log((l+e)./(l-e)); 

Ao  =  ((l-(e.A2))./(e.A3)).*(log((l+e)./(l-e))  -  2*e); 
kl  =  Ao./(2  -  Ao); 
k2  =  Bo./(2  -  Bo); 

kb  =  ((e.A4).*(Bo  -  Ao))./((2  -  e  A2).*(2*e  A2  -  (2  -  e  A2).*(Bo  -  Ao))); 

%  Normal  force  coefficients 

%  Drag  coefficient  determination 

for  i  =  1:102  -  k 
if  ln/d(i)  <  1 

Cdo(i)  =  cln(l)*(ln/d(i))A4  +  cln(2)*(ln/d(i))A3  +  cln(3)*(ln/d(i))A2  +  ... 
cln(4)*(ln/d(i))  +  cln(5); 

else 

Cdo(i)  =  0.29; 
end 

if  lb(i)/d(i)  <  1 

Cdo(i)  =  clb(  1  )*(lb(i)/d(i))A4  +  clb(2)*(lb(i)/d(i))A3  +  ... 
clb(3)*(lb(i)/d(i))A2  +  clb(4)*(lb(i)/d(i))  +  clb(5); 

else 

Cdo(i)  =  0.29; 
end 
end 

lv  =  0.905  *1B; 

r  =  (d./2).*(l  -  (lv  -  If) ./lb); 
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Sv  =  pi.*r.A2; 

Cla  =  (2*(k2  -  kl).*Sv)./Sb; 

Yvprime  =  (-Sb./lB.A2).*(Cla  +  Cdo); 

Yvp  =  [Yvp  Yvprime  -(ones(l,k-l))]; 

Cma  =  (2*(k2-kl)./(Sb*lB)).*((-pi.*d.A2./12).*(3*xm  -  In)  +  ... 
((pi.*d.A2)./(12*lb.A2)).*  ... 

(-3*xm.*lv.A2  +  6*lv.*xm.*lb  +  2*lv.A3  -  3*lf.*lv  A2  -  3*lb.*lv  A2  +  . 
6*lv.*xm.*lf  -  6*lf.*xm.*lb  +  3*lb.*lf  A2  -  3*xm.*lf.A2  +  If  A3)); 

Nvprime  =  (Sb.*Cma)./lB.A2; 

Nvp  =  [Nvp  Nvprime  -(ones(l,(k-l)))]; 

%  Rotary  force  coefficients 

Stb  =  r.*(lB  -  lv); 

Cmq  =  Cma.*((l  -  xm./lB)  A2  -  V*(lcb  -  xm)./(Stb.*lB  A2))./... 

((1  -  xm./lB)  -  (V./(Stb.*lB))); 

Nrprime  =  -(Sb.*Cmq)./lB.A2; 

Nrp  =  [Nrp  Nrprime  -(ones(l,(k-l)))]; 
clear  Nrprime 

Clq  =  Cla.*(l-  xm./lB); 

Yrprime  =  (-Sb.*Clq)./lB  A2; 

Yrp  =  [Yrp  Yrprime  -(ones(l,(k-l)))]; 

%  Acceleration  coefficients 

Zwdotprime  =  (2*k2*V)./(lB  A3); 

Yvdotprime  =  -Zwdotprime; 

Yvdp  =  [Yvdp  Yvdotprime  -(ones(l,(k-l)))]; 
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for  I  =  1 : 102-k; 


Iydf_ln(I)  =  quad8('vdn',0,ln); 

Iydf_lm(I)=  quad8('vdm',ln,lf(I)); 

IydfJb(I)  =  quad8('vdb',lf(I),lB); 
end 

Iydf  =  Iydf_ln  +  Iydf_lm  +  Iydfjb; 

Nrdotprime  =  (-2.*kb.*Iydf)./(rho.*lB.A5); 

Nrdp  =  [Nrdp  Nrdotprime  -(ones(l,(k-l)))]; 

N  =  isnan(Nrdp); 

I  =  find(N  >  0); 
for  t  =  1  :length(I) 

Nrdp(I(t))  =  -1.0; 
end 

clear  d  If  In  lm  lb  xm  I 

clear  Iydf_ln  Iydf_lm  Iydf_lb  Iydf 

clear  global 

IB  =  4.9941; 

rho  =  62.4/32.174; 

end 

Nvdp  =  -Yvdp; 

Nvdp  =  reshape(Nvdp,101,101); 

Yvp_d  =  reshape(Yvp,101,101); 

Nvp_d  =  reshape(Nvp,101,101); 

Yrp_d  =  reshape(Yrp,101,101); 

Nrp_d  =  reshape(Nrp,101,101); 

Yvdp_d  =  reshape(Yvdp,101,101); 

Nvdp_d  =  reshape(Nvdp,101,101); 

Nrdp_d  =  reshape(Nrdp,101,101); 

save  Yvp_d  Nvp_d  Yrp_d  Nrp_d  Yvdp_d  Nvdp_d  Nrdp_d 
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%  Format  for  presentation  (all  analysis  complete) 

Nvp  =  reshape(Nvp,101,101); 
figure(l) 

plot(R,Nvp(l :  101,6), 'c',  R,Nvp(l:101,ll),'r+',  R,Nvp(l:101,16),’b*',  ... 

R,Nvp(l:101,21),'go’) 
xlabel('Mid-body  Fraction') 
ylabel('Nvprime') 

text(0. 1,-0. 0132, 'Nose  Fraction  (%)') 
legend('  5','10','15','20') 
axis([0  1  -0.0145  -0.013]) 
hold  on 

plot(R,Nvp(l:  101,1  l),'r:’,  R,Nvp(l:101,16),'b:',  R,Nvp(l:101,21),'g:') 

hold  off 

pause 

print  -depsc2  nvp 

Nvp  =  reshape(Nvp,  1,10201); 
k  =  find(Nvp  >  -0.01); 

Nvp(k)  =  -0.02*(ones(length(k),l)); 

Nvp  =  reshape(Nvp,  101,1 01); 

k  =  find(Nrp  <  -0.01); 

Nrp(k)  =  ones(length(k),l); 

Nrp  =  reshape(Nrp,101,101); 

figure(2) 

plot(R,Nrp(l:101,6),'c',  R,Nrp(l:101,ll),'r+',  R,Nrp(l:101,16),'b*',  ... 

R,Nrp(l:101,21),'go') 
xlabel('Mid-body  Fraction') 
ylabel('Nrprime') 

text(0. 2, 0.0025, 'Nose  Fraction  (%)') 
legend('  5','10','15',’20') 
axis([0  1  -0.0025  0.003]) 
hold  on 

plot(R,Nrp(l:  101,1  l),’r:',  R,Nrp(l:101,16),’b:',  R,Nrp(l:101,21),'g:’) 

hold  off 

pause 

print  -depsc2  nrp 
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Nrp  =  reshape(Nrp,  1,10201); 
k  =  find(Nrp  >  0.05); 

Nrp(k)  =  -0.02*(ones(length(k),l)); 

Nrp  =  reshape(Nrp,101,101); 

figure(3) 

plot(R,Nrdp(l:101,6),'c',  R,Nrdp(l:  101,1  l),’r+',  R,Nrdp(l:101,16),'b*',  ... 

R,Nrdp(  1:101,21  ),'go') 
xlabel('Mid-body  Fraction’) 
ylabel('Nrdotprime') 
text(0. 6, -0.00015, 'Nose  Fraction  (%)') 
legend('  5’, ’10’, ’15', ’20’) 
axis([0  1  -0.0012  0]) 
hold  on 

plot(R, Nrdp(l:  101,1  l),'r:',  R,Nrdp(l:101,16),'b:’,  R,Nrdp(l:101,21),'g:') 

hold  off 

pause 

print  -depsc2  nrdp 

k  =  find(Nrdp  <  -0.02); 

Nrdp(k)  =  -0.02*(ones(length(k),l)); 

Nrdp  =  reshape(Nrdp,101,101); 

k  =  find(Yvp<-0.02); 

Yvp(k)  =  -0.02*(ones(length(k),l)); 

Yvp  =  reshape(Yvp,101,101); 

k  =  find(Yrp  <  -0.02); 

Yrp(k)  =  -0.02*(ones(length(k),l)); 

Yrp  =  reshape(Yrp,101,101); 

k  =  find(Yvdp  <  -0.02); 

Yvdp(k)  =  -0.02*(ones(length(k),l)); 

Yvdp  =  reshape(Yvdp,101,101); 

figure(4) 

plot(R,Yvp(l :  101 ,6),'c',  R,Yvp(l:  101,1  l),’r+’,  R,Yvp(l:101,16),'b*',  ... 

R,Yvp(l:101,21),'go') 
xlabel('Mid-body  Fraction') 
ylabel('Yvprime') 

text(0. 15, -0.0075, 'Nose  Fraction  (%)') 
legend('  5','10’,'15','20’) 
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axis([0  1  -0.01  -0.004]) 
hold  on 

plot(R, Yvp(  1 : 1 0 1 , 1 1  ),'r R,Yvp(l:101,16),’b:',  R,Yvp(l:101,21),'g:’) 

hold  off 

pause 

print  -depsc2  yvp 
figure(5) 

plot(R,Yrp(l:101,6),'c’,  R,Yrp(l:101,ll),’r+',  R,Yrp(l:101,16),'b*',  ... 

R,  Yrp(  1:101,21  ),'go’) 
xlabel('Mid-body  Fraction') 
ylabel('Yrprime') 

text(0. 15, -0.0025, 'Nose  Fraction  (%)') 
legend('  5', '10', '15', '20') 
axis([0  1  -0.005  0]) 
hold  on 

plot(R,Yrp(l:  101,1  l),'r:',  R,Yrp(l:101,16),'b:',  R,Yrp(l:101,21),'g:') 

hold  off 

pause 

print  -depsc2  yrp 
figure(6) 

plot(R,Yvdp(l :  101 ,6),'c',  R,Yvdp(l:101,l l),'r+',  R,Yvdp(l:101,16),'b*',  ... 

R,Yvdp(l:101,21),'go') 
xlabel('Mid-body  Fraction') 
ylabel('Yvdotprime') 
text(0. 5, -0.0146, 'Nose  Fraction  (%)') 
legend(’  5','10','15’,'20') 
axis([0  1  -0.0154  -0.0145]) 
hold  on 

plot(R,Yvdp(l:101,ll),'r:',  R,Yvdp(l:101,16),'b:',  R,Yvdp(l:101,21),’g:’) 

hold  off 

pause 

print  -depsc2  yvdp 
figure(7) 

plot(R,Nvdp(l :  101 ,6),'c',  R,Nvdp(l :  101,1  l),'r+’,  R,Nvdp(l:101,16),... 

'b*',  R,Nvdp(l:  101,21), 'go') 
xlabel('Mid-body  Fraction') 
ylabel(’Nvdotprime') 
text(0. 4,0. 015, 'Nose  Fraction  (%)') 
legend('  5',’10’,'15’,'20') 
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axis([0  1  0.0145  0.0154]) 
hold  on 

plot(R,Nvdp(l:101,ll),'r:',  R,Nvdp(l:10U6),'b:',  R,Nvdp(l:101,21),'g 

hold  off 

pause 

print  -depsc2  nvdp 
%  Range  change  for  mesh  graphics 
R  =  (0:0.02:1); 
figure(8) 

S  =  reshape(S,101,101); 
mesh(R,R,S(  1:2: 101, 1:2: 101)), grid 
xlabel('Nose  Fraction') 
ylabel('Mid-body  Fraction') 
zlabel('Slenderness  Ratio') 
view(60,30) 
print  -depsc2  slender 

figure(9) 

mesh(R,R,Yvp(  1:2: 10 1,1:2: 101)),  grid 

xlabel('Nose  Fraction') 

ylabel('Mid-body  Fraction') 

zlabel('Yvprime') 

view(60,30) 

print  -depsc2  yvpm 

figure(10) 

mesh(R,R,Nvp(l:2: 101,1:2:101)),  grid 

xlabel('Nose  Fraction') 

ylabel('Mid-body  Fraction') 

zlabel('Nvprime') 

view(40,40) 

print  -depsc2  nvpm 

figure(l  1 ) 

mesh(R,R,Yrp(l:2:101, 1:2:101)),  grid 
xlabel(’Nose  Fraction') 
ylabel('Mid-body  Fraction') 
zlabel('Yrprime') 
view(60,30) 
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print  -depsc2  yrpm 


figure(12) 

mesh(R,R,Nrp(l:2: 101,1:2:101)),  grid 

xlabel('Nose  Fraction') 

ylabel('Mid-body  Fraction') 

zlabel('Nrprime') 

view(40,30) 

print  -depsc2  nrpm 

figure(13) 

mesh(R,R,Yvdp(l:2:101, 1:2:101)),  grid 

xlabel('Nose  Fraction') 

ylabel('Mid-body  Fraction') 

zlabel('Yvdotprime') 

view(60,30) 

print  -depsc2  yvdpm 

figure(14) 

mesh(R,R,Nvdp(l:2:101,l:2:101)),  grid 

xlabel('Nose  Fraction') 

ylabel('Mid-body  Fraction') 

zlabel('Nvdotprime') 

view(20,20) 

print  -depsc2  nvdpm 

figure(15) 

mesh(R,R,Nrdp(l:2:101,l:2:101)),  grid 

xlabel(’Nose  Fraction') 

ylabel('Mid-body  Fraction') 

zlabel('Nrdotprime') 

view(60,50) 

print  -depsc2  nrdpm 

%  Clipping  of  data  for  graphical  presentation 

Yvp_s  =  reshape(Yvp,  1,1 0201); 
k  =  find(  Y  vp_s<-0 . 01); 

Yvp_s(k)  =  -0.01*(ones(length(k),l)); 

Yvp_s  =  reshape(Yvp_s,101,101); 

Nvp_s  =  reshape(Nvp,  1,10201); 
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k  =  find(Nvp_s  >  -0.013); 

Nvp_s(k)  =  -0.0145*(ones(length(k),l»; 
k  =  find(Nvp_s  <  -0.0145); 

Nvp_s(k)  =  -0.0145*(ones(length(k),l)); 
Nvp_s  =  reshape(Nvp_s,101,101); 

Yrp_s  =  reshape(Yrp, 1,10201); 
k  =  find(Yrp_s  <  -0.005); 

Yrp_s(k)  =  -0.005*(ones(length(k),l)); 
Yrp_s  =  reshape(Yrp_s,101,101); 

Nrp_s  =  reshape(Nrp,  1,1 0201); 
k  =  find(Nrp_s  <  -0.005); 

Nrp_s(k)  =  -0.005*(ones(length(k),l)); 
Nrp_s  =  reshape(Nrp_s,101,101); 

Yvdp_s  =  reshape(Yvdp, 1,10201); 
k  =  find(Yvdp_s  <  -0.0154); 

Yvdp_s(k)  =  -0.0154*(ones(length(k),l)); 
Yvdp_s  =  reshape(Yvdp_s,101,101); 

Nvdp_s  =  -Yvdp_s; 

Nrdp_s  =  reshape(Nrdp,  1,10201); 
k  =  find(Nrdp_s  <  -0.00125); 

Nrdp_s(k)  =  -0.00125*(ones(length(k),l)); 
Nrdp_s  =  reshape(Nrdp_s,  10 1,101); 

figure(16) 

mesh(R,R,Yvp_s(l:2:101,l:2:101)),  grid 

xlabel('Nose  Fraction') 

ylabel('Mid-body  Fraction') 

zlabel('Yvprime') 

view(60,30) 

print  -depsc2  yvps 

figure(17) 

mesh(R,R,Nvp_s(l  :2: 101 , 1:2:101)),  grid 
xlabel('Nose  Fraction') 
ylabel('Mid-body  Fraction') 
zlabel('Nvprime') 
view(40,40) 
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print  -depsc2  nvps 


figure(18) 

mesh(R,R,Yrp_s(l:2:101,l:2:101)),  grid 

xlabel('Nose  Fraction') 

ylabel('Mid-body  Fraction') 

zlabel('Yrprime') 

view(60,30) 

print  -depsc2  yrps 

figure(19) 

mesh(R,R,Nrp_s(  1:2: 10 1,1:2: 101)),  grid 

xlabel('Nose  Fraction') 

ylabel('Mid-body  Fraction') 

zlabel('Nrprime’) 

view(40,40) 

print  -depsc2  nrps 

figure(20) 

mesh(R,R,Yvdp_s(l:2:101, 1:2:101)),  grid 

xlabel('Nose  Fraction') 

ylabel('Mid-body  Fraction') 

zlabel('Yvdotprime') 

view(60,30) 

print  -depsc2  yvdps 

figure(21) 

mesh(R,R,Nvdp_s(l  :2: 101 ,1 :2: 101)),  grid 

xlabel('Nose  Fraction') 

ylabel('Mid-body  Fraction') 

zlabel('Nvdotprime') 

view(40,30) 

print  -depsc2  nvdps 

figure(22) 

mesh(R,R,Nrdp_s(  1:2:101,1:2:101)),  grid 

xlabel('Nose  Fraction') 

ylabel('Mid-body  Fraction') 

zlabel('Nrdotprime') 

axis([0  1  0  1  -0.00125  -0.0005]) 

view(60,50) 

print  -depsc2  nrdps 
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B.  VDN.M  PROGRAM 


%LCDR  Eric  Holmes 
%Thesis 

%Iydf_ln  Functional  Integration 

%  d  =  diameter 
%  IB  =  total  body  length 
%  In  =  nose  length 
%  lm  =  middle  length 
%  lb  =  base  length 

%  xm  =  geometric  middle  of  the  body 
%  rho  =  density  of  water 

function  II  =  vdn(x) 

global  d  IB  If  In  lm  lb  xm  rho  I 

11  =  ((-rho*pi.*d(I).A2.*x.*(xm  -  x).A2./(4.*lnA2)).*(x  -  2.*ln)); 

C.  VDM.M  PROGRAM 

%LCDR  Eric  Holmes 
%Thesis 

%Iydf_lm  Functional  Integration 

%  d  =  diameter 
%  IB  =  total  body  length 
%  In  =  nose  length 
%  lm  =  middle  length 
%  lb  =  base  length 

%  xm  =  geometric  middle  of  the  body 
%  rho  =  density  of  water 

function  12  =  vdm(x) 

global  d  IB  If  In  lm  lb  xm  rho  I 

12  =  (rho*pi.*d(I).A2/4).*(xm  -  x).A2; 
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D.  YDB.M  PROGRAM 


%LCDR  Eric  Holmes 
%Thesis 

%Iydf_lb  Functional  Integration 

%  d  =  diameter 
%  IB  =  total  body  length 
%  In  =  nose  length 
%  Im  =  middle  length 
%  lb  =  base  length 

%  xm  =  geometric  middle  of  the  body 
%  rho  =  density  of  water 
%  If  =  fore  body  (lb  +  lm) 

function  13  =  vdb(x) 

global  d  IB  If  In  lm  lb  xm  rho  I 

13  =  (pi*rho.*d(I).A2.*(xm  -  x)  A2).*(lb(I)  -  x  +  lf(I))  A2./(4*lb(I)  A2); 
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APPENDIX  B.  SAMPLE  MATLAB  ROUTINE  FOR  DETERMINING 
FUNCTIONAL  COEFFICIENTS  FROM  A  SURFACE  PROFILE 


A.  SURF.M  PROGRAM 

%  LCDR  Eric  P.  Holmes 
%  Surface  Equation  Solver 

n  =  2; 

x  =  [];  y  =  []; 

load  hcdata 

Yvp  =  Yvp_d(6:26, 41:61); 
Yvp  =  Yvp(:); 

%  Range 

xr  =  (5:25); 
for  i  =  1:21 
x  =  [x;xr]; 
end 

x  =  0.01  *x; 
x  =  x(:); 


yr  =  (40:60); 
for  i  =  1:21 
y  =  [y  yr]; 

end 

y  =  0.01*y; 

y  =  y(0; 


%  Program  operation 

n  =  n+1; 
k  =  1; 

A  =  zeros(size(x)); 
for  i  =  n:-l:l, 
for  j  =  1  :i 

A(:,k)  =  ((x.A(i-j)).*(y.A(]-l))); 
k  =  k+1; 
end 
end 
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p  =  (A\Yvp).' 

%  Error  calculation 


x  =  reshape(x,21,21); 
y  =  reshape(y,21,21); 

Yvp  =  reshape(Yvp,21,21); 

for  i  =  1:21 
for  j  =  1:21 

Yvpcal(ij)  =  p(l)*x(i,j)A2  +  p(2)*x(i,j)*y(i,j)  +  p(3)*y(i,j)A2  +  . 
p(4)*x(i,j)  +  p(5)*y(i,j)  +  p(6); 
end 
end 

for  i  =  1:21 
for  j  =  1:21 

perr(ij)  =  (Yvpcal(ij)  -  Yvp(i,j))/Yvp(i,j); 

end 

end 

perr  =  100*perr; 
xp  =  x(l,:); 

yp  =  y(:,i)’; 

figure(l) 

mesh(xp,yp,Yvp),  grid 
xlabel('Nose  Fraction') 
ylabel('Mid-body  Fraction') 
zlabel('Yvprime') 
view(60,30) 
print  -depsc2  yvp 

figure(2) 

mesh(xp,yp,perr),  grid 
xlabel('Nose  Fraction') 
ylabel('Mid-body  Fraction') 
zlabel('Yvprime  Percent  Error') 
axis([.05  .25  .40  .60  -20  20]) 
view(60,30) 
print  -depsc2  yvperr 
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